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Abstract 

We analyze the possibilities of triplet pairing in Sr2Ru04 based upon an idea of 
interlayer coupling. We have considered two models differing by the effective interactions. 
In one model the quasi-particle spectra have horizontal line nodes on all three Fermi 
surface sheets, while in the other the spectra have line or point nodes on the a and 
/3 sheets and no nodes on the 7 sheet. Both models reproduce the experimental heat 
capacity and penetration depth results, but the calculated specific heat is sightly closer 
to experiment in the second solution with nodes only on the a and j3 sheets. 

1. Introduction. 

Strontium Ruthenate (S^RuCU) is widely believed to be a spin triplet superconductor^ 
121 Hill], however the theoretical model and particular pairing mechanism are still hotly de- 
bated (51 IB1 171 151 l9l ITDl 111! 1121 113| . Its lattice structure resembles the layered structure of 
the cuprate La2- : rBa :E Cu04 but, instead of high temperature d-wave superconductivity, the 
superconducting phase of strontium ruthenate (T c ~ 1.5K) appears to be p-wave or f-wave in 
nature. The strongest evidence for this comes from the 17 O NMR Knight shift data[H| and 
neutron scattering experiments^] which indicate that the in-plane Pauli spin susceptibility 
is constant below T c . These experiments could be naturally explained with a triplet pairing 
state d(k) = e z (k x ± ik y ) in exact analogy with the ABM phase of superfluid 3 He. The 
/iSR experiments 16, also indicate a spontaneous breaking of time reversal symmetry at T c , 
which would be consistent with such a chiral ABM type state. On the other hand several 
experiments|17l 1181 119| indicate that the gap function must have lines of nodes on the Fermi 
surface, unlike the simple ABM state which is nodeless on the three cylindrical sheets, a, /3 
and 7, of the Fermi surface of Sr2Ru04. 

In this paper we address the question of the possible location of these line nodes on the 
Fermi surface. We concentrate on the case of horizontal lines of nodes, because vertical nodes 
(for example in f-wave pairing states) appear to be inconsistent with the absence of angular 
dependence of the thermal conductivity in an a-b plane magnetic field 19 20 . The presence 
of horizontal lines of nodes, as originally suggested by Hasegawa, Machida and Ohmi 5 , 
cannot be explained in any 2-d theoretical model but requires a 3-d model with at least 
some component of the pairing interaction acting between planes. A number of such models 
assuming interlayer coupling have been proposed |d1 171 151 l9l 1101 [TTI 1121 
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The specific question which we address here is whether there are horizontal line nodes on all 
three Fermi surface sheets, as proposed recently by Koikegami, Yoshida and Yanagisawa| 1 1 1 . 
or whether the nodes are only on the a and j3 sheets, as proposed by Zhitomirsky and 
Rice[Sj and in our earlier interlayer coupling model|51 1101 112| . To answer this question we 
will analyze, an effective week coupling model where the attractive interactions can appear 
between electrons on nearest neighbour and next nearest neighbour lattice sites (symbolised 
by 2 and 3 lines in Fig. 1). We will contrast the predictions of our original interlayer coupling 
model[21 11UI I12| with ones chosen to reproduce the gap structure proposed by Koikegami, 
Yoshida and Yanagisawa|l 1 1 . 




Figure 1: Body-centred tetragonal lattice of Ruthenium atoms in the Sr2Ru04 structure. 
The full lines show possible interactions between electrons occupying the single site '1' and 
its in-plane nearest neighbour sites '2' and intra-plane neighbours '3'. 
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2. The interlayer coupling model 

To describe superconductivity in S^RuCU, we start from the following simple multi-orbital 
attractive Hubbard Hamiltonian, 



ijvnm' ,<r 

2 ^ ] U mm i {ifyflirnofljm' a > . 



(1) 




Figure 3: Temperature dependence of order parameters A^ 1 , A C c X (a) and A^ 1 and 
A^ c x (b) found for different choice of interactions. Fig. 3a corresponds to the results for 
only U i (ij) < (Eq. 2 - scenario 1) while Fig. 3b to the case where U , (ij) < for 

m, m! = a,b (Eq. 3 - scenario 2) and J/J,c(*j) < 0. 



Here i and j label the sites of a body centred tetragonal lattice (as shown in Fig. 1), and m 
and m! refer to the three Ruthenium t2 g orbitals (as shown in Fig. 2). In the following the 
orbitals will be denoted a = xz, b — yz and c = xy. The hopping integrals t mm i (ij) and site 
energies e m were fitted to reproduce the experimentally determined three-dimensional Fermi 
surface |21l I22| . We found that the set t mm i given in |12|gave a good account of the Fermi 
surface data. 

Since the actual physical mechanism of pairing in S^RuCU is unknown, we will adopt a 
phenomenological approach and treat the effective Hubbard interaction constants U™ m i (ij) 
as free parameters. The idea is to test different "scenarios" for the interaction constants 
against the available experimental results. Good agreement between experiment and theory 
is likely to only occur when the interaction parameters U^ m i (ij) lead to a gap function on 
the Fermi surface d(k) which is similar to that actually present in the material. In particular, 
both the jump in specific heat at T c , and the linear dependence of C(T)/T (C(T)/T ~ T 
for line nodes) near to T — are sensitive to the gap d(k) over the whole Fermi surface. 
Scenarios in which the line nodes occur on all Fermi surface sheets or only on a would 
be expected to lead to different temperature dependencies of C(T), which can therefore be 
distinguished by comparison to the experiments. 

In Fig. 1 we depict schematically the possible interactions U™ m , (ij) between electrons 
on the same Ru-atom (shown by circle T') and the neighbouring Ru-atoms (shown by lines 
denoted '2', '3' for in-plane and out of plane interactions, respectively). We assume that 
the dominant interaction between electrons on the same site i (T') is the strong Coulomb 
repulsion. The effective interactions between nearest neighbours are assumed to be attrac- 
tive. They could arise either from spin-fluctuation mediated exchange^, or from interlayer 
Coulomb scattering|ll|. 

In this paper we will compare two different model sets of interaction constants. The 
first is motivated by the inter-plane Coulomb scattering model of Koikegami, Yoshida and 
Yanagisawa 11 .. It assumes Coulomb repulsion between electrons on the same and the nearest 
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Figure 4: Minimum energy quasiparticle eigenvalues on the a (a) , /3 (b) and 7 (c,d) Fermi 
surface sheets. Fig. 4c (scenario 1) corresponds to the results for only £/ x m /fo) < while 

Fig. 4d (scenario 2) to the case where U 1 - < for m,m' = a,b and C/J jC (ij) < 0. 

Fig. 4a and b are the same for both cases. 



neighbour in plane lattice sites ('1' and '2' in Fig. 1) and attraction between Ru planes ('3' 
in Fig. 1). The corresponding Hubbard interaction parameters are 

/ \ / u± u ± \ 

U 11 = and U x = I U± U 1 - , (2) 

V J \ u ,J - J 

expressed as matrices in orbital space, m, m' . In this model there are two attractive Hubbard 
parameters, one acting only between c orbitals and the other acting equally between a and b 
orbitals. In k-space there is little hybridization between the a-b orbitals and the c orbitals, 
and so these inter-plane interactions mainly corresponding to interactions within the 7 band 
{U' ± ) and within the a and (3 Fermi surface sheets (U±). We shall call this case 'scenario 1' 
in the rest of this paper. 

We wish to compare the predictions of the above set of model parameters with those of our 
previous inter-plane coupling model|§l 1101 112| . which was motivated by the different spatial 
orientations of the xz, yz, xy orbitals as shown in Fig. 2. Given that the Ru d-xy orbital (c) 
has a mainly 2-d character we assumed that interactions between c orbitals are mainly in- 
plane. On the other hand the Ru d-xz (a) and d-yz (b) orbitals are oriented perpendicular to 
the RUO2 plane, and so we assumed that the interactions between electrons in these orbitals 
are mainly iner-plane. This simple reasoning leads to the following two parameter Hubbard 
model with 

/ \ / u ± u ± \ 

U 11 = and U x = t/- 1 U 1 - . (3) 

\00t/ [l / \ / 




We shall refer to this as 'scenario 2' below. 

In our earlier papers |5] I I1JI I12| we showed that this simple model gives a good overall 
account of the temperature dependent heat capacity C(T), in-plane superfluid density n s (T), 
and thermal conductivity. In |12| we showed that the predictions of the model are fairly 
robust against the addition of extra interaction parameters or disorder. 

In equations 2 and 3 we have set to zero any interaction terms which are either assumed to 
be small, or those which may be repulsive. We have checked that all the zero values appearing 
in the above interaction matrices Ij" and U 1 " (Eqs. 2,3) can be changed into small positive 
values representing repulsions without any change to the solution, and so we decided to use 
(Eqs. 2,3) the minimal set leading to pairing. 

For above choices of interactions within the negative U extended Hubbard model (Eq. 1), 
we solved the Bogolubov-de Gennes equations: 

[ E»-B m , m ,(ij) A™' m ,(ij) \(u» m „, 
{ A*^(ij) E» + H mm , ) \ v» m , a , 

together with the self-consistency condition 

C'Km't')); Xmm'(ij) = u ima v jm'a'(l ~ 2 /(-E"')), (5) 

which follow from Eq. (1) on making the usual BCS-like mean field approximation ,22]. Here 
f{E v ) is the Fermi function, j3 = 1/fcsT, ks is Boltzmann constant and v enumerates the 
solutions of Eq. 0] 

Assuming an p-wave pairing state of the form d(k) ~ e z , on each Fermi surface sheet, 
then we need only consider the gap parameters A , (k) at each point in the Brillouin zone. 




A" 

1 Ti 




Figure 6: Calculated superfuid density, function of temperature, T, compared to 

the experimental data of Bonalde et aZ.|18l (gray squares). The line '1' shows the results 
for only r(ij) < while line '2' corresponds to the case where U < for 

m, m! = a,b and £/J,e(*i) < 0. 



Dropping the spin indices for clarity, the general structure of pairing parameter is of the 
general form 

A mm - (k) = A w sin k x + A v mm , sin k y 

. \ z • k z c k x ky . t . k x . ky . k z c 
+A mm / sm cos — cos — + A , sm — sm — sm 

mm 2 2 2 mm 2 2 2 

/ a X • ky . A y ky ky \ k Z C , . 

+ A„, m / sm — cos — + A , cos — sm — cos (6) 

I mm mm' o o / o \ ' 



for m, m' — a, b and c. In the present calculations we neglected the possibilities of p z pairing 
(A^ m ,) or f-wave pairing A^ m , , for reasons which are discussed further in |51ll2|. 



3. Line Nodes and Specific Heat 



From the two different interaction models given by Eqs. 2 and 3 we numerically find the 
corresponding solutions of the gap equation Eqs. 4, 5. In the case of scenario '1', Eq. 2, the 
gap parameters have the general form, 

Aw (k) = ( A x mm , sin y cos y 1 + Af nm , cos y sin y j c °s ^y > (7) 

for m, m' = a, b, c. While for scenario '2' (Eq. 3) the gap parameters are 
A cc (k) = A* c sin k x + A^ c sin k y 

A rom / (k) = I A w sm -y cos y + A^ m , cos — sm y I cos — (8) 
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for m,m' — a or b, respectively. Clearly scenario '1' has a order parameter for which all of 
the gap parameters A mm /(k) vanish in the planes k z — ±7r/c2. Therefore all three Fermi 
surface sheets should have horizontal line nodes. On the other hand, in scenario '2' only the 
a and b components of A mm / (k) vanish, implying that the 7 sheet is nodeless. 

The temperature dependence of the order parameters in each scenario are shown in Fig. 
3. In our analysis we have fitted the interaction parameters (U" = — 0.494i, U ± = — 0.590t, 
U' ± = — 0.312t) to obtain a single critical temperature T c « 1.5K for all order parameter 
components (Fig. 3). Note, the differences between U ± = — 0.590t and U' ± = — 0.312£ arises 
as an effect of the difference in partial density of states for the bands a, f3 and 7. In case of 7 
the Fermi surface is close to a Van Hove singularity |11II24| . and so U' ± can be smaller while 
still obtaining the same T c . On the other hand, the larger value of U can be explained by 
the out of plane spatial orientation of the a and b Ru orbitals, compared to the in plane Ru 
d-xy orientation of c orbital 

The angular dependence of the eigenvalues Ek itP , i — x, y on the corresponding p = a, f3, 
7 Fermi surface sheets, plotted along the z axis, are shown in Fig. 4. In case of the inter- 
plane attraction only scenario '1' the gap has line nodes on all three Fermi surface sheets. In 
contrast, in scenario '2' the gap is nodeless on the 7 sheet, as can be seen in Fig 4(d). 

Now the key question is whether experiment can distinguish between these two gap sce- 
narios, namely line nodes on all Fermi surface sheets compared to just nodes on a ft only. 
We calculated the specific heat for those solutions via the following relation: 

k,p 

where / is the Fermi function. 

The results are presented in Fig. 5 and compared to the experimental values of Nishizaki 
et aZ. [T7|. 

One can see in Fig. 5 that the slope of C(T)/T near to T = is slightly higher for scenario 
T' compared to scenario '2', consistent with the "extra" line node on the 7 Fermi surface 
sheet. However the change in slope is quite small, and so one can say that either model is 
consistent with the low temperature experimental data. 

On the other hand, in Fig. 5 it is clear that the second solution '2' works slightly better for 
the jump of specific heat at critical temperature T c . Solution '1' has a smaller Fermi surface 
average of |A(k)| 2 , and this leads to the slightly smaller jump in specific heat at T c . The close 
similarity of curves T' and '2' can be understood if one notices that the quite different gap 
symmetry on the 7 sheet can lead to rather similar results after integration over the Fermi 
surface (Eq. 10). 

As a further test of the presence of horizontal line nodes on all Fermi surface sheets, we 
have also calculated the in-plane superfluid density n s (T) using 




where X(T) denotes the temperature dependent penetration depth, Vi = v x or v y is the in 
plane band velocity at k, e is the electron charge, no = 4-7T x 10~ 7 is the magnetic constant, 
£kp is the electron band energy and 

MT) = A>)_ 

n s (0) A 2 (T)- 1 ' 

Both secenarios '1' and '2' give similar results of n s and agree with experimental results 
by Bonalde et aZ. |18| (gray squares). However, one can note that the slightly different slope 
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in small temperature regions gives a slight advantage to solution '2' which mimics the exper- 
imental data a little better. 

4. Conclusions. 

We have tested two different gap models for strontium ruthenate, which are consistent 
with two physically different pairing mechanisms. In scenario '1', we assumed that all in- 
plane interactions are repulsive, and that only the out of plane interactions lead to pairing. 
This is motivated by the Coulomb scattering pairing mechanism of Koikegami, Yoshida and 
Yanagisawa'll . In scenario '2' we assumed attractive in-plane interactions for Ru d — xy 
orbitals (7 band) and attractive inter-plane interactions of the a and b orbitals (a and f3 
bands). Surprisingly we found that the predicted specific heat is very similar in both models 
(Fig. 5), even though one has a horizontal line node on all three Fermi surface sheets while 
the other has a nodeless 7 sheet. Similarly the temperature dependent superfuid density is 
closer to experiment in both scenarios. Of these two models the nodeless 7 sheet appears to 
be slightly close to the experiment, but the actual differences are small. 

In these calculations we have not attempted to include the interband proximity effect 
proposed by Zhitomirsky and Rice|5j. In |12| we showed that this corresponds in real-space 
to the addition of three-site Hubbard interaction parameters, or assisted hopping, to the 
Hubbard Hamiltonian. It would be quite possible to combine this interband proximity effect 
interaction with either of the two pairing scenarios which we have considered here. However 
the close similarity of the specific heat and superfluid density in the two models, shown in 
Figs. 5, 6 strongly suggests that the effect of the proximity coupling terms would also be 
very similar in either pairing model. 
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